Limits on the WIMP-nucleon scattering cross-section from neutrino telescopes 

G. Wikstrom^B and J. Edsjo^Q 

1 The Oskar Klein Centre for Cosmoparticle Physics, Dept. of Physics, 
AlbaNova, Stockholm University, SE- 10691 Stockholm, Sweden 
(Dated: March 18, 2009) 

Neutrino-telescopes like Super-Kamiokande and IceCube have started to explore the neutrino 
fluxes from WIMP annihilations in the Sun. The non-observation of a signal can put constraints on 
the WIMP properties. We here focus on the neutrino signal from WIMP annihilation in the Sun and 
show that under reasonable assumptions, the non-observation of a signal from IceCube puts a much 
tighter constraint on the spin-dependent WIMP-proton scattering cross-section than current direct 
detection experiments like COUPP and KIMS. For the spin-independent scattering cross-section, 
the limits from IceCube and current direct detection experiments like XENONIO and CDMS place 
similar constraints. We here go through the assumptions being made and the uncertainties that arise 
in converting from limits on the muon flux from the Sun to limits on the WIMP-proton cross-section, 
and present our results as easy to use conversion factors. 



I. INTRODUCTION 



One of the favourite candidates for the dark matter in the Universe is Weakly Interacting Massive Particles 
(WIMPs), of which neutralinos that arise in supersymmetric extensions of the standard model is a popular example 
(see e.g. Ref. |l[ for a review). If WIMPs are the dark matter in the Universe, they will be abundant in the Milky Way 
halo where they can produce annihilation radiation (gamma rays, neutrinos, positrons, antiprotons, antideutrons, . . . ) 
that can be searched for with e.g. the Pamela Q and Fermi Q satellites. The WIMPs can also scatter in heavy objects 
like the Sun [ij], where they can lose enough energy to be captured gravitationally and eventually sink to the core. 
Once in the core, they can annihilate and produce e.g. neutrinos that can be searched for with neutrino telescopes 
like Super-Kamiokande ||, IceCube @ and Antares Q. Another route to search for the WIMPs is direct detection 
experiments where one tries to observe WIMP scatterings on a detector target. For both the neutrino signal from the 
Sun and the direct detection experiments, the strength of the signal depends crucially on either the spin-independent 
and/or the spin-dependent WIMP-nucleon scattering cross-section. Currently, XENONIO [8[ and CDMS @ place the 
most stringent direct limits on the spin-independent scattering cross-section, whereas KIMS [l(j and COUPP [ll| 
put the most stringent direct limits on the spin-dependent scattering cross-section. However, as the neutrino signal 
from the Sun as searched for by neutrino telescopes also depends strongly on these cross-sections, we can use these 
searches to put limits on the same scattering cross-sections. This is of course more involved than for direct detection 
experiments as we need to consider the capture in the Sun, the subsequent annihilation, production of neutrinos, 
propagation, interaction and oscillation of neutrinos, and finally scattering near the neutrino telescope to produce 
observable muons (/i + and Most of these steps can be done in a fairly straightforward approach under reasonable 
assumptions, and the limit derived on the spin-dependent scattering cross-section is highly competitive with those 
from direct detection experiments. A similar analysis was done in Ref. 12], but those results need updating, mainly 
because neutrino oscillations change the neutrino fluxes, but also because better modelling of hadronization/decay 
of quark jets affect the results. We will here perform a more detailed study, including full three flavour neutrino 
oscillations. For this, we will use the WimpSim [r| WIMP annihilation Monte Carlo as built into DarkSUSY [l4[ 
5.0.3, and derive conversion factors to go from limits on the muon flux from neutrino telescopes to limits on the 
spin-independent and spin-dependent WIMP-proton scatte ring cross-sections. At the end, we will apply these on 
the recent results from Super-Kamiokande [f| and IceCube [l5f . The correlation between signals in direct detection 
experiments and neutrino telescopes was also investigated in Ref. fl6|. in an analysis based on the results in [l2j , but 
adding annual modulation to the picture. In another recent paper [17| . the correlation between neutrino telescopes 
and direct detection experiments was revisited for light WIMP dark matter (below 20 GeV). We will here be more 
general than in that paper and focus on a wide range of WIMP masses. 
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II. WIMP CAPTURE AND NEUTRINOS 



As there are many assumptions needed to calculate the neutrino flux from the scattering cross-section, we will 
do this in two ways. First, we will focus on a 'standard' calculation, with reasonable assumptions usually made in 
the literature. We will then go through some of these assumptions and perform a conservative calculation where we 
instead change our assumptions to see how low fluxes we can reasonably get for given scattering cross-sections. 



A. Standard calculation 



The number of neutralinos in the Sun, N, is described by the differential equation 

^ = C c - C A N 2 - C E N, (1) 

where the three constants describe capture (Cc), annihilation (Ca), and evaporation (Ce). For the masses of interest 
here, the evaporation is small and can be neglected fl8l |. Note that Cc depends on the scattering cross-sections on 
the elements in the Sun, whereas Ca depends on the annihilation cross-section. The neutralino annihilation rate Ta 
is then 

Y A = \C A N 2 = iCctanh 2 (t/r) (2) 



r=\I^Cc~C~ A , (3) 

where the present rate is found for t = t Q ~ 4.5 • 10 9 years. When t Q /t ^ 1 annihilation and capture are in 
equilibrium, — 0, and 



dt 

Ta = \c c - (4) 

Thus, in equilibrium, Y a only depends on the capture rate Cc, i.e. Ta only depends on the scattering cross-sections, 
and not on the annihilation cross-section. From Ta it is then straightforward (but tedious) to calculate the resulting 
muon flux at a neutrino telescope for a given WIMP annihilation channel. 

For the capture (Cc) and annihilation (Ca) we follow the treatment in Ref. as implemented in DarkSUSY [Til ]. 
In short, we have to integrate the WIMP scatterings over the volume of the Sun and the velocity of the WIMPs in 
the galactic halo and will here review the most essential steps in this calculation. 

For the velocity of the WIMPs, we will here assume that it follows a Maxwellian distribution, 

where u is the velocity of the WIMP (outside the potential well of the Sun), vq — 220 km/s is the velocity of the 
Sun relative to the halo, Vd = 270 km/s is the WIMP velocity dispersion, and n x is the WIMP number density. We 
assume that the local WIMP density is 0.3 GeV/cm 3 . At a given interaction point in the Sun, the WIMP velocity is 
given by w = y/u 2 + v 2 , where v is the escape velocity at that point. 

The WIMP capture rate per unit shell volume from element i in the Sun is [l^l 

dC . ru ma 

— — = / du wll Vt i{w), (6) 

dV Jq u 

where Q v i (w) is the capture probability per unit time for element i. In Eq. (O we integrate up to the velocity 
Umax = 2«y / /I/(/i — 1) at which the WIMPs scatter to the escape velocity v. Here /x = m x /rrii, for the WIMP mass 
m x and the mass m, of element i. 

For scattering on heavier elements than Hydrogen, we also need to take decoherence into account by introducing a 
form factor suppression. Following standard lore, we assume a Helm-Gould exponential nuclear form factor (l9l | for 
the momentum transfer q on element i 

iF^^expi-AE/E?) ; AE = ^— (7) 
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where 



3fi 2 

2m v -R?' v ; 



0.91(m 2 /GeV) 1/3 + 0.3 
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Qv,i{w) can then be calculated analytically and we arrive at the expression [lj 

(^ + 1) 2 



tov,i(w) = a X i ■ m — — — E t 



6XP 1 ~^Ef ) ~ 6XP \T^+W mX ^Ef 



2n u 2 + v 2 



(9) 



Here a X i is the WIMP scattering cross-section on element i, and rij is the number density of element i in the Sun. 
To get the total capture rate we integrate over the radius of the Sun, _Rq, and finally arrive at [T^ | 

Co- r^dr^^, (10) 







where the solar composition is taken from Ref. [2Cj|], where we take the solar model BS2005-OP as our default choice. 
For elements heavier than Oxygen, the solar model BS2005-OP only gives the total abundance and we use the relative 
abundances as given in Ref. [2l| for these. In total, we include the 16 most important elements up to Nickel. 

7 X A, 



We can calculate the spin-independent (SI) cross section on a nucleus with mass number A, cr^, as pj 



a si * IA * {m x m A f (m x + m p y 

a xA — G A 7 1 Y2~t Y2~ ^ > 

A [m x + niA) \m x ra v Y 

where a SI is the SI WIMP-proton scattering cross-section. 

For the spin-dependent (SD) WIMP-proton cross-section a SD , Hydrogen is the most important element in the Sun, 
and we can neglect the heavier elements in this case. 

The muon flux and at the detector is calculated from the WIMP annihilation rate in the Sun as 

*" = Sj5? f ^ f ^ f° dA r < p ^K> ^ da ^;; K) e p(* o e ( 12 ) 

47ri^ j E th j E th Jo j Efi ai/^ '-r' *y ati u 

where Dq is the distance from the center of the Sun, n the target number density, E 1 * the detector's energy threshold, 
A the muon range, P(E fll £J„, A) the probability for a muon of energy E^ to have a final energy E^ after a path-length 
A in the detector material, da v {E v , E^/dE^ is the differential neutrino cross-section for production of a muon with 
energy E^ from a neutrino of energy E u , P(fi,i) is the probability that a produced neutrino of flavour i oscillates 

to flavour /z in the detector, Bf is the branching ratio for annihilation channel /, and dNj jdE v is the differential 
number of neutrinos of flavour i produced per annihilation in channel /. We perform this calculation with WimpSim 
[l3| including interactions, neutrino oscillations and regeneration of neutrinos from tau decay. We have performed 
these simulations with the neutrino oscillation parameters [11] (see [l3| for details) Q\2 = 33.2°, 6i 3 — 0°, 623 — 45.0°, 
5 = 0, Am^ = 8.1 x 10~ 5 eV 2 and lAm^l = 2.2 x 10 _3 eV 2 but our results are not very sensitive to the exact choice 
of these parameters. The results of these WimpSim simulations are included in DarkSUSY [l4[ as easy to use tables 
and interpolation routines. For our purpose here, we use a muon energy threshold of 1 GeV as this is how most 
neutrino telescopes report their results. We will also use the muon fluxes without any angular cut on the muons. 

We then have all the necessary ingredients to calculate the flux of muons (or neutrinos) in a neutrino telescope from 
any given WIMP scattering and annihilation cross-sections. However, for the neutrino telescope limits to be really 
competitive, we need to have a situation where we have equilibrium between annihilation and capture as discussed 
above (i.e. we need to have T cq < 4.5 • 10 9 years). If we do not have equilibrium, neutrino telescope limits will 
still constrain WIMP models, but we lose the direct correlation between the neutrino-induced muon fluxes and the 
scattering cross-sections. To investigate this further, we need to look at specific WIMP models and we will here focus 
on the neutralino that arises as a natural WIMP dark matter candidate in supersymmetric extensions of the standard 
model for particle physics, like the Minimal Supersymmetric Standard Model (MSSM). We have made extensive scans 
of the MSSM model parameters both in mSUGRA models and low-energy phenomenological MSSM-7 and MSSM-9 
models (in these models, we specify the supersymmetric parameters at the electroweak scale instead of at the GUT 
scale as in mSUGRA, see [23[ for details). In Fig. Q]we show the neutrino- induced muon flux <5> p for models passing 
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current accelerator constraints and with a relic density in the cosmologically preferred range (generously chosen), 
together with recent limits from IceCube 15], and Super-Kamiokande [a], and the full IceCube estimated sensitivity 
[151 ] . In this figure we indicate for which models we do have equilibrium between capture and annihilation and for 
which we do not. For all models with high fluxes that are currently probed and will be probed in the near future 
by neutrino telescopes, equilibrium is fulfilled and we will from now on assume that we do have equilibrium. For 
hypothetical scenarios where this assumption would be invalid, it is straightforward to rescale our results with the 
tanh 2 (i/r oq ) factor given in Eq. ^ above. We also note that for another popular class of models, Kaluza-Klein 
dark matter in Universal Extra Dimensions (UED) [24| . we also always have equilibrium or near equilibrium between 
capture and annihilation. Naively, we do expect to have equilibrium for models that can be probed with currently 
planned experiments, as the annihilation cross-section that enters r in Eq. © is constrained by the relic density 
constraint to be of the order of (crv) ann ~ 3 • 1CP 27 cm 3 s _1 . The above examples for MSSM and UED confirm that 
this is the case for these specific models. 
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FIG. 1: Muon flux from the Sun as a function of neutralino mass, for mSUGRA and MSSM models where neutralinos are 
(dots) or are not (crosses) in equilibrium in the Sun. Also shown are limits from IceCube [l5[, and Super-Kamiokande [fj, and 
the full IceCube estimated sensitivity fl~5j ] . The muon flux is given above a threshold of 1 GeV. 



B. Conservative calculation and uncertainties 



The calculation above includes 'standard' assumptions both from particle physics and astrophysics. However, we 
want to use the neutrino telescope results to put a limit on the scattering cross-section. Hence, we need to investigate 
processes that can affect our results. In particular we want to investigate the processes that can reduce the neutrino 
fluxes since a reduction would mean that our limits on the scattering cross-sections would be increase (i.e. they would 
have been too optimistic). 

First of all, there are some uncertainties in the solar model, especially regarding the abundances of heavy elements. 
In a conservative setup, we instead focus on the solar model BS2005-AGS,OP [2(| which includes new lower estimates 
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of the abundances of the heavier elements (for the actual relative abundance we use the estimates in [21| for heavier 
elements than Oxygen). This will not affect capture via spin-dependent scattering as this occurs on Hydrogen, 
whose abundance is well-known. It will however reduce capture via spin-independent scattering, especially for heavier 
WIMPs (heavier than a few hundred GeV), where capture on the heavier elements is important. 

Another effect that can reduce the capture rate is gravitational effects from the planets, especially Jupiter, during 
the capture process [2^, [26[. In Ref. [IE, it was estimated that for low scattering cross-sections, the life time of 
WIMP orbits after the first scattering in the Sun is quite long (i.e. it takes a long time before they scatter again and 
eventually sink to the core). If the WIMP orbit reaches out to Jupiter, there is a large probability that Jupiter will 
affect the orbit and disturb it so that it no longer passes through the Sun and eventually throws the WIMP out of the 
solar system back into the Milky Way halo. To get a conservative estimate of this effect, we have modified the capture 
expressions above so that WIMPs that reach out to Jupiter after the first scatter in the Sun will be considered as not 
captured. This reduces capture, especially for heavy WIMPs and especially for scattering on light elements in the 
Sun (i.e. it is more pronounced for spin-dependent scattering). For 1 (10) TeV WIMPs the reduction is about 13% 
(87%) for spin-dependent capture and about 1% (12%) for spin-independent capture. 

The form factors introduce another source of uncertainty in the spin-independent capture, and the usually assumed 
Helm-Gould form factors are rather approximate (but convenient to use). In Ref. [2jJ, the effect of different form 
factors on direct detection rates was investigated and the conclusion was that the total scattering rate could change 
by up to about 20% with more accurate form factors. For our conservative estimate, we will here use the Helm-Lewin- 
Smith form factors [27, 28], which for mass number A are given by 

R = yj(c? + (7/3)7r 2 a 2 - 5s 2 ) (13) 
s = 0.9 • 10~ 15 m, a ~ 0.52 • 10~ 15 m, c ~ [1.23.4 1/3 - 0.60] • 10~ 15 m. 

Another source of uncertainty is of course the local dark matter density and velocity distribution. For our conser- 
vative estimate, we keep the local halo density fixed to po = 0.3 GeV/cm 3 as for the standard calculation, as it is easy 
to rescale our results for other densities. The local WIMP density also enters in the same way for direct detection 
experiments, so it does not affect the conversion factors below. A bigger concern is the velocity distribution as di- 
rect detection experiments are more sensitive to high velocity WIMPs (as they produce larger energy nuclear recoils), 
whereas the capture process in the Sun is more efficient for low velocity WIMPs. N-body simulations indicate (see e.g. 
[29l |) that dark matter halos probably are both anisotropic and deviates from the simple Maxwellian form. However, 
dynamical constraints [3(| indicate that they cannot change by too much and still be consistent with observations. It 
seems reasonable to assume though, that the velocity distribution is not that different from a Maxwellian distribution, 
that it would reduce the capture in the Sun by more than a factor of two. For our conservative estimate, we have 
though not included the effects of changing the velocity distribution, as there is yet no consensus on how different it 
could be and still be compatible with observations. There is also a recent proposal [3l[ that there might be a dark 
matter disk in the Milky Way in addition to the regular dark matter halo. This disk could have a local density of 
(0.25-1.5) po and would corotate with the stars in the galaxy. The velocity distribution of this population would 
therefore be peaked at much lower velocities than the regular halo, which would cause very little change in direct 
detection experiments (except if the energy thresholds could be reduced), but would increase the capture rate in the 
Sun by maybe one order of magnitude [321 ] . The existence of this dark matter disk is still debated and as it increases 
the capture rate instead of decreasing it, we do not include it in our conservative estimate. We just want to point out 
that if such a disk exists, the neutrino telescopes could be even more sensitive than the standard calculation gives. 

III. RESULTS 

We are now ready to calculate our results and will first focus on the conversion rates from muon flux limits to 
scattering cross-section limits. We will then compare with earlier calculations and at the end we will apply our 
conversion factors on the current best limits from neutrino telescopes and compare with direct detection experiments. 

A. Conversion factors 

If we assume that one annihilation channel dominates (i.e. assuming Bt — 1 for some annihilation channel / in 
Eq. (fT2"|) ) we can relate the observed muon flux to the annihilation rate 



(14) 
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using the functions 77^ (m x ), shown for four (extreme) channels in Fig. [2] These values are the same in the standard 
and conservative calculations. Note that for a specific WIMP model, we typically know the branching fractions 
Bf, but for general WIMPs, we do not have this knowledge and in extreme cases the uncertainty coming from this 
lack of knowledge could be about an order of magnitude as seen in the Fig. [5J However, for most WIMP models, 
the uncertainty is usually much smaller than this. For supersymmctric neutralinos, the branching fraction to hard 
channels is typically > 10% for m x > mw, so the sum of 77's is typically rather close (within a factor of a few) to the 
values for a hard channel like W + W~ . For Kaluza-Klein dark matter in models of Universal Extra Dimensions, the 
branching fractions are rather well-known and with a large hard contribution. Hence, for most realistic models, ry is 
closer to the hard channels like W + W~ than the soft ones like bb. To get a really conservative limit, we could use the 
soft bb channel as an extreme case, but we have here chosen to present results for these different channels separately 
Further on, to get a conservative limit on either the spin-independent or spin-dependent scattering cross-section, 
we assume that either one of these dominates capture. We can then through Eqs. (fTT)|) and get the WIMP-proton 
scattering cross-sections a SI or a SD as a function of Ta 

<t si = X SI (m x )r A (15) 
a SD = \ SD (m x )T A , 

for functions X SI (m x ) and X SD (m x ). We can then finally relate the cross-section to the muon flux, as 

Kf(m*)$£ (16) 

K fV x )$/. 

The conversion factors Kf(m x ) and K^ D (m x ) between cross-section and muon flux are shown for four annihilation 
channels in Fig. [3] for the standard and conservative calculations. For the SI case, we note that the main effect of the 
conservative calculation is to increase k si due to the different form factor. At higher masses, the lower abundances of 
heavy elements and the effect of Jupiter increases slightly further. For the SD case, we see essentially no differences 
between the two calculations at low masses as the only significant difference between the calculations in this case is 
the reduction of capture due to Jupiter, that starts increasing k sd at high masses. 

The conversion factors in Fig. [3] can be used to convert a muon flux limit from a neutrino telescope to a limit on 
the SI or SD WIMP-nucleon scattering cross-section. We have also developed a web tool [33| that can be used to 
get accurate values for these conversion factors for any WIMP mass between 10 GeV and 10 TeV. The web tool 
also contains more annihilation channels than those in the figure and also allows to convert from fluxes with different 
muon energy thresholds and angular cuts than those assumed here. The web tool is also available as a downloadable 
Fortran program to be used with DarkSUSY [l4| • 
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FIG. 2: The conversion factor - J£ fr^pr — as a function of ra x , for annihilation channels / = W + W~ , t + t~ , it, and bb. This 
conversion is identical in the standard and conservative calculations. 
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annihilation channels / = W W , r r , tt, and bb. Lower lines for each channel are from the standard calculation and 
upper lines are from the conservative calculation. For the muon flux we have used a threshold of 1 GeV. 



B. Comparison with earlier calculations 

We can now compare our results with earlier estimates. In Kamionkowski et al. the event rates in neutrino 

telescopes were compared with the expected event rates in a hypothetical Hydrogen direct detection detector (and 
were given in units of m kg ). They assumed a set of different branching fractions Bf for different mass ranges 
and gav e their result as the spin-dependent scattering rate. Using the same expression for the direct rate, given in 



Ref. [34(, we have calculated the equivalent conversion factors Kj{m x ) in terms of the cross-section instead of the 
direct detection event rate. In Fig. [4] we show these conversion factors from Ref. 1131 (dotted line) and compare 
with our results with our standard calculation (solid line). As the calculation in Ref. [12j does not include neutrino 
oscillations, for comparison reasons we also show our results without including neutrino oscillations (dashed line). 

In this plot, the branching fractions for the black lines are assumed to be / = bb for m x < mw — 80 GeV and 
/ = W + W~ or for m x > 80 GeV. The grey lines are instead for / = t+t~ for m x < m t ~ 175 GeV and / = tt 

for m x > m t ~ 175 GeV. In this figure, we see some discrepancies with the calculation by Kamionkowski et al.. 

First of all, our conversion factors for bb are higher than their results, which is mostly due to a different treatment 
of hadronization and decay of the b quarks. We have used Pythia [35| for hadronization/decay, whereas [l2| have used 
analytical expressions that e.g. do not include final state radiation. Also, the treatment of b meson interactions in the 
Sun is slightly different in our treatment where we have used updated interaction cross-sections (see [l|| for details). 
Our results for bb agree well with those in Ref. [36l ]. and as they are softer than in 12] our conversion factors are 
higher. For W + W~ and tt, our results agree fairly well, but for t + t~, we see a large discrepancy between the grey 
solid and dashed curves (below m t ). This discrepancy comes from neutrino oscillations not being included in the [12| 
treatment and our results with neutrino oscillations turned off (grey dashed curve) , agree fairly well with the results in 
[l2l ]. The reason neutrino oscillations have such a big effect on the t + t~ channel is that we get significantly more tau 
neutrinos than other neutrino flavours from the r decays and oscillations will efficiently transform these to detectable 
muon neutrinos at the detector. Hence, the neutrino flux is significantly enhanced due to neutrino oscillations, which 
will result in a lower conversion factor in Fig. 2] 

Other smaller differences come from us using an updated solar model and the use of the full expressions for solar 
capture (including an integration over the solar radius) . 

We have also compared with the calculation in Ref. [T3|, where they focus on WIMP masses up to 20 GeV. They 
have re-evaluated the Super-Kamiokande limits for this low-mass case and used an updated method of the one in [l2j . 
Their limits on the cross-section are lower than those we have here, where we use the Super-Kamiokande published 
results (see Fig. [5]). In [ItJ they analyze the Super-Kamiokande data up to 30° from the Sun, where there happens to 
be a downwards fluctuation in the number of recorded events compared to the background expectation, which results 
in a better limit. They have also used approximate expressions for WIMP capture in the Sun and similar analytic 
expressions for the neutrino yield from bb as in [l2l ]. These effects taken together can explain the difference between 
our results. 
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FIG. 5: Limits on spin-dependent WIMP-proton cross-section as a function of neutralino mass, for Super-K calculated from the 
muon flux limit [jj, and for IceCube-22 Limits from CDMS 0, COUPP KIMS [3, and the full IceCube estimated 
sensitivity [l5| are shown for comparison. Also shown are mSUGRA and MSSM models where neutralinos are (dots) or are 
not (crosses) in equilibrium in the Sun 



C. Limits on scattering cross-sections 

Using our conversion factors in Fig. [3] (we here use the standard calculation), we can now convert limits on the muon 
flux from neutrino telescopes to limits on the scattering cross-section. As expected, the limits on the spin-independent 
scattering cross-section are are similar to, but not competitive with, current direct detection limits and we will instead 
focus on the spin-dependent scattering cross-sections. In Fig. [5] we show the limits on the spin-dependent scattering 
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cross-section for our standard calculation using the recent muon flux limits from the IceCube 22-string detector [la 
and the full IceCube estimated sensitivity [l5j . For comparison, we also show the limits from Super-Kamiokande [f|. 
For Super-Kamiokande, we convert their muon flux limit from a threshold of 1.6 GeV to limits on the flux with a 1 
GeV threshold. They have also applied an angular cut in the range of 1.5° to 19° [37| for different masses, and we 
also correct for this effect to get a limit on the full muon flux without angular cuts. The limits are further given for a 
mixture of 80% bb, 10% cc, 10% t + t~ , and we have here chosen to calculate the limit for bb. In practice, the choice of 
annihilation channel mixtures in the Super-Kamiokande analysis mainly affects the choice of angular cuts, and as we 
have corrected for that, it is a good approximation to treat the limits as limits on the flux for any annihilation channel 
(the limits on the cross-section will depend on the annihilation channel though through the channel dependence in 
K S f D ). M 

For comparison, in Fig. [5l we also show the direct detection limits on the spin-dependent scattering cross-section 
from CDMS COUPP [ll|, and KIMS [lj3]. As can be seen, the neutrino telescope limits are highly competitive, 
even if we take the annihilation channel uncertainties into account. Even if we would switch to the conservative 
calculation, the limits from neutrino telescopes would be far better than current direct detection experiments. 



IV. CONCLUSIONS 



We have calculated the conversion factors between muon fluxes in neutrino telescopes and direct detection scattering 
cross-sections and applied these to the current neutrino telescope limits from IceCube and Super-Kamiokande. The 
limits on the spin-dependent scattering cross section are shown to be highly competitive with current direct detection 
experiments. 

The main difference between our result and the previous result of Kamionkowski et al. [T2| comes from including 
neutrino oscillation effects in the Sun. For / = bb, there is another difference, mainly coming from how b hadroniza- 
tion/decay and 6-meson interactions are treated. In addition there is a smaller difference coming from a more exact 
calculation of the capture in the Sun in our treatment. 

Finally, using these new conversion factors, experiments like IceCube, Super-Kamiokande and Antares can easily 
convert their limits on the muon flux to limits on the scattering cross sections (as was recently done by IceCube in 

El). 
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For the Super-Kamiokande limits, there is a discrepancy between our calculated limits on a SD and the limits Super- 
Kamiokande themselves have calculated [f|] on a SD (by using the results in [i"^). Their limits are about a factor of 5 higher 
and this difference cannot be explained with our (small) differences with the results in fl2l|. It probably stems from an 
inconsistent use of different direct detection event rate calculations when using the results in 1^]. We have been in contact 
with the authors of p] about this possible error. 



